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1 Introduction and Scope. 

o\ : 

Q\ • In the solid state physics of crystals, an important part is played by the 

band structure of the electronic states. This band structure arises from 
the representation of the periodicity in the electronic state space. Powerful 
J ■ computational methods were developed for the calculation of band structures, 

among them the linear muffin-tin (LMTO) method |l^J in the atomic sphere 
approximation (ASA) In the physics of quasicrystals, it is believed that 
>- ! the electronic system plays an important part ||i~9| . Here one is lacking the 

periodic symmetry. To still use the powerful methods of band computations, 
one must replace the quasicrystal by a periodic approximant. The question 
arises how such approximant computations approach a quasiperiodic limit. In 
a recent calculation it is shown by a supercell analysis that an approximant 
computation may lead to artefacts in the electronic density of states (DOS). 
More detailed local properties of the electrons appear as derivatives of the 
electronic charge density in Mossbauer studies on quasicrystals |T^| . 
In the present work we wish to analyse and compare electronic states in 
periodic and quasiperiodic potentials in a way free of approximations. In 
this way we can hope to address the similarities and differences of these 
systems from first principles. 

To a first and very important approximation, the many-electron states in or- 
dered crystalline or quasicrystalline solids are constructed from one-electron 
states, compare for example Ashcroft and Mermin fl2|, after appropriate an- 
tisymmetrization. The one-electron Hamiltonian then contains periodic and 
quasiperiodic potentials respectively. For a general account of band theory 
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we refer to Blount ||]. Analytic properties and the approach due to Wan- 
nier [2(J are discussed by Kohn ||. For general electronic structures beyond 
periodicity, a local view was advocated by Heine [[/]] who proposed to throw 
out K-space. 

In what follows we wish to elaborate and to compare electronic systems and 
their local structure in terms of the following concepts: 

(i) Local view of if-space in periodic potentials: 

A Hamiltonian with an infinite periodic potential has the discrete symmetry 
group A of its lattice. The electron states are characterized by the irreducible 
representations of A. These are labelled by the continuous set of Bloch vectors 
K from the Brillouin zone, that is, the first unit cell of the reciprocal lattice 
A^. The eigenstates for fixed energy E(K) are the Bloch states belonging to 
the bands. 



A local formulation of K-space arises as follows: Wigner and Seitz pi 



m 



1934 introduced the cellular method: The one-electron Schrodinger equation 
with fixed energy E is to be solved exclusively on the finite unit cell of the 
lattice with the local boundary conditions that the solution, after propagation 
over a full primitive period say aj, j = 1,2,3, picks up a pure phase factor. 
This phase factor when written in the form fj := exp(iKjaj), j = 1, 2, 3 de- 
termines the Bloch labels. Two solutions with K — > —K are degenerate. The 
energy E(K) appears in bands, with Bloch labels ranging over the Brillouin 
zone, and in gaps where there are no states with Bloch type boundary con- 
ditions. This approach allows an extension to finite and to infinite systems. 
Within bands, the solutions can be matched and propagated as Bloch states 
over all the cells of the crystal. 

(ii) Local view of i^-space in quasiperiodic potentials: 

For infinite almost periodic systems in ID some general results with many 
references are discussed in |J. For another detailed study of discrete systems 
with many references we quote Siito fl~8|| . For our purposes we consider an 
infinite quasiperiodic potential built from a few basic cells as in a tiling model 
for quasicrystals. For each cell one can apply the local analysis of Wigner 
and Seitz with local Bloch type boundary conditions. Quasiperiodicity now 
requires a matching of these local boundary conditions between the basic 
cells. It is not enough that, at fixed energy E, each cell separately admits 
Bloch type boundary conditions. We must also check if a patch of cells admits 
these boundary conditions. Again one can study the quasiperiodic extension 
of the system, following for example a geometric inflation of the tiling. 

(iii) Local view, positive and negative energy, motive clusters: 

Imagine the atomic nuclei of a unit cell in a periodic crystal to form a free 
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motive cluster. This motive cluster is not unique since the unit cell of a 
crystal admits translations as well as transformation of shape, for example 
in going from the primitive to the Wigner-Seitz cell. 

For a fixed choice of motive clusters we should not relax the positions of the 
nuclei to gain energy, as would be done in real molecules, since in the crystal 
the positions are controlled by all atomic neighbours. 

Example: Take for example a linear crystal with two atoms A, B, in two 
spacings a ^ b with the sequence AaBbAaBbAa... and period a + b. Choose 
two different motive clusters by the replacements AaBb — > AaB and BbAa — ■» 
BbA which do not affect the bound states. We get in general two different 
bound states, with the binding energy depending on the relative distances 
a and b respectively. In a linear molecule we would relax the distance and 
end up in both cases at the same binding energy. This is not allowed in 
the crystal! In the linear crystal the two different motive clusters generate 
identical periodic structures and bandsl 

Consider then the one-electron bound states for the chosen motive cluster. 
It is true that, in order to treat it like a molecule, we should take along the 
nuclei and their Coulomb repulsion. However since we keep the nuclei at 
fixed positions we may neclect these interactions which of course influence 
the binding and the total energy. Independent of these interaction terms, the 
different boundary conditions of the electron states at positive or negative 
energy will enter into the computation of the many-electron states. 
Since the Hamiltonian for fixed nuclear positions is the same as in the Wigner- 
Seitz method, the only difference is in the boundary conditions which for a 
bound state require the solutions to decay exponentially. Clearly now we 
must distinguish between the positive and negative energy region. The former 
admits scattering states, the latter bound states. 

The infinite systems at positive energy should here be considered as an ap- 
proximation to a finite macroscopic systems which still admits a scattering 
experiment. We could also consider half-infinite systems which allow at least 
for backscattering. For both finite and infinite systems, the distinction be- 
tween positive and negative energy states is crucial. In considering finite 
patches of periodic or quasiperiodic systems one therefore should pay atten- 
tion to the sign of the energy. 

(iv)Bound, Bloch and scattering states in periodic potentials: 

We shall term the negative energy region the exact tight binding scenario in 
what follows. As a partial differential equation of second order, the stationary 
Schrodinger equation has two fundamental solutions on a single cell. Suppose 
that the energy of a bound state lies within the range of energy for a band. 
Then it follows from the Schrodinger equation that the bound state solution 
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restricted to the unit cell must be an exact linear combination of the local 
Bloch states at the bound state energy. In the positive energy region we 
can still obtain states with local Bloch boundary conditions. Now we can 
inquire how the band structure of these positive energy states is related to 
one-electron scattering states. 

Example: The relation of bound and Bloch states in crystals is often consid- 
ered in approximate rather than exact computations, as for example in the 
tight binding approximation. Here one starts exclusively from bound atomic 
1-orbitals of the isolated atoms, obtains the Bloch states from their tunneling 
to next neighbours, and labels the bands by the 1-orbitals. If first-order de- 
generate perturbation theory applies, the energies in the bands result from 
the splitting of the single orbital whose degeneracy is proportional to the 
number of atoms. In this approximation it is clear that the energy of the 
local atomic orbital is within the corresponding band. In many cases the 
bound atomic orbitals of the atoms in the cell are separated in energy. Oth- 
erwise one can consider hybrid or motive cluster orbitals on the unit cell as 
the origin of band labels. The discrete version of this approximation leads 
to an eigenvalue problem and therefore looses the information on the sign of 
the energy. 

The general arguments given above are valid independent of any such ap- 
proximation scheme. They show that in the exact tight binding scenario 
there must exists an exact local relation between bound and Bloch states. 

(v)Bound, Bloch and scattering states in quasiperiodic potentials: 

For a quasiperiodic potential based on a tiling, we again wish to distinguish 
the negative and positive energy regions. One can consider at negative en- 
ergy the bound states for finite patches which form part of a quasiperiodic 
structure and compare with states obeying local Bloch type boundary con- 
ditions on the same patch. We refer to such an energy range as a band germ 
and to the local Bloch states as Bloch germs. If the bound state energy 
admits Bloch type boundary conditions, it must be possible to express the 
bound states by pairs of local Bloch states (germs). The bound states may 
be related to local clusters. 

In the scattering from quasiperiodic potentials at positive energy, the aim is 
again to compute the scattering matrix. For an exact computation one may 
inquire about the implication of inflation symmetry on the scattering matrix. 

In what follows we wish to demonstrate the validity of these concepts in a 
way free of approximations. In dimension three we cannot yet implement 
exact examples for these concepts. The LMTO method may allow to exam- 
ine some of our points. We choose here a comparative and specific study 
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restricted to continuous periodic and quasiperiodic potentials in dimension 
one. On them we wish to demonstrate the concepts given above by exact 
computations. To do so we pay the price of choosing simple delta potentials 
even of the same strength. We explore the electron states in terms of local 
boundary conditions, in the spirit of the Wigner-Seitz approach, applied to fi- 
nite strings (patches) which form part of periodic or quasiperiodic structures. 
We shall always start from finite patches and then look into their recursive 
extension. The solutions with Bloch boundary conditions on a finite string 
are systematically compared with bound state solutions on the same string. 
The main tool of our analysis is the continuous transfer matrix which is 
well defined for piecewise constant potentials, including delta potentials as a 
limiting case. This matrix propagates by matrix multiplication a fundamen- 
tal system of two solutions. We shall stress in what follows the polynomial 
dependence of the matrix elements on the strength parameter of the delta 
potential. This will allow us to order and analyze the matrix systems. 



2 Finite periodic strings at negative energy. 

2.1 Preview: An energy gauge for crystals. 

As an introduction we present some material from fl2||. We rephrase the well- 



known periodic case fll5|| . Consider first a finite string S with the transfer 
matrix M, described in more detail in section 2.2. We define an energy 
gauge f as a function with value / = on an energy interval such that 
|tr(M)| < 1 and / = 1 otherwise, compare Fig. 1. We call this a band 
germ. The condition on the trace assures that inside the band germ M has 
two complex conjugate eigenvalues of absolute value 1. Repeat the string 
n times to produce a new transfer matrix M n . Since the existence of band 
germs is related to the eigenvalue problem, the range of energies for which 
^\tr(M n )\ < 1 is independent ofn, the band germs stay the same on the new 
string. The Bloch germs propagate through the string and pick up the same 
phase factors respectively after each transmission. 

Consider next the bound states of the string M n . In a very tight binding 
case we claim 

1 Prop: The bound states of the string M n may be grouped into sets of n 
states, where the energies of each such set corresponds to a part / = of the 
energy gauge / of the string M. 

Proof. For very tight binding it suffices to use first-order degenerate per- 
turbation theory, applied to the bound states of the single atoms: To this 
order, the bound states of the string are the eigenvalues of a matrix whose 
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off-diagonal entries are the weak atom-atom cross terms. By standard matrix 
theory, the maximum level splitting will increase with n. But band theory 
tells us that in the limit n — > oo all energies stay inside the band, hence 
inside the initial band germ. Thus the energy of all these bound states for 
finite n must stay within the energy gauge of the (single) band germ. 
We now have the following situation in the finite string M n : In a band germ 
from the string M , there is for each energy value a pair of Bloch germs which 
can carry charge current. There are now n discrete bound states with the 
energy gauge as in the initial string. 

f(S) S f{S 2 ) S 2 f{S A ) S A 



Fig.l: Periodic strings: The string S to the left has one attractive 5- well with 
a single bound state, followed by a tunnel. The vertical bar to its left shows 
the energy interval for the band germ. This interval is the energy gauge f(S) 
comprising the bound state. The string S 2 in the middle has two attractive 
(^-potentials and two bound states. The energy gauge f(S 2 ) is unchanged 
but comprises two bound states. The same energy gauge /(S* 4 ) for the string 
S 4 to the right comprises four bound states. 



A schematic view of the periodic scheme is given in Fig.l. Now we can extend 
the analysis to n — > oo: 

2 Prop: For an infinite periodic repetition of a fixed string, the energy gauge 
stays the same as for the initial string. The band germ generates a band. 
Within each band there is an infinite set of (pairs of) Bloch states whose 
energies Ek < fill up the original band germ. 



2.2 Bloch and bound states in a single band. 

We shall use the notation of Consider the ID Schrodinger stationary 
equation on the line x G R, 



h^d^_ 

2m dx 2 



+ V{x) 



(x) = E<f>(x). 



(1) 
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with V(x) a piecewise constant potential. Denote two fundamental solutions 
and their derivatives by <fii(x), 4>i(x), 2 (x),0' 2 (x). Choose the initial values 
0i (0) = 1, 4>i(0) = 0, 02(0) = 0, 02 (0) = 1 which implies a Wronskian equal 
to 1. The 2x2 standard transfer matrix is defined as 



M(x) = 
M(0) = 1. 



02 (^) 



X 0n(X 



(2) 



The transfer matrix obeys the first-order system of equations which is equiv- 
alent to the Schrodinger equation. A second interpretation of the transfer 
matrix, which we shall adopt in what follows, results because the matrix 
M(x) can also be shown to propagate the fundamental system. 
The Schrodinger equation for piecewise constant potentials on a finite string 



can now be solved as follows p2| : We first determine the transfer matrix for 
simple building blocks of finite size at the fixed energy. Then we propagate 
the solutions over the full string by matrix multiplication. This matrix mul- 
tiplication guarantees the continuity of the solutions and of their derivatives. 
Specific boundary conditions can finally be met by determining the appropri- 
ate linear combination of the two fundamental solutions. For real matrices 
with Wronski determinant 1 the transfer matrices belong to the matrix group 
SX(2, R). More details on this group and the isomorphic group £77(1, 1) can 
be found in |flO| . One should be careful in applying the group concepts to 
the transfer matrices, as can be seen from the following transformation prop- 
erties. 

It will prove convenient to pass to a new fundamental systems of solutions 
by choosing different initial values. With respect to the matrix formed by 
the solutions, this transformation is achieved by right multiplication with a 
constant matrix, 

M(x) -> M(x)C (3) 

The new transfer matrix which propagates the new system M(x)C is no longer 
equal to M(x)C, instead it is given by C~ 1 M(x)C . This propagating transfer 
matrix by itself does not admit an interpretation in terms of fundamental 
solutions and their derivatives! 

We now pass to a new system of fundamental solutions with different initial 
conditions by right multiplication with the matrix 



R 




1 1 

— K K 



(4) 



The system corresponding to the matrix M(x)R has the initial data R at 

2m 



x = 0. At negative energy E = — we define a tunnel as a string with 
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vanishing potential. In a tunnel the transfer matrix takes the form 




exp(— kx) 
— Kexp(— kx) 



exp(/txj 



(5) 



The transfer matrix which propagates this new system at negative energy is 
given by 



M(x) := R- 1 M{x)R 

kx) 

exp(nx) 



(6) 



exp( 




For a bound state we require that an exponentially increasing function at neg- 
ative energy on the left-hand tunnel, after passing an intermediate transfer 
matrix M, produces an exponentially decreasing function on the right-hand 
side. In the new basis, this property and the form eq. |5] require that the 
intermediate transfer matrix obeys M22 = 0. 

3 Prop: In the system M of transfer matrices for finite strings, the condition 
for a bound state requires that the second diagonal element vanishes. 
A tunnel of length b, followed by an attractive delta-potential of strength u 
at negative energy E = — we term the string S. We define 5 := u/k 
and Ai := exp(/3), (3 := Kb to obtain the transfer matrix of S as 



-¥ (1-5*) 

Aifl- 



Ar 1 






Ai 



(7) 



-Ar 1 ^ 



¥) 



Here the first transfer matrix describes the delta-potential and the second 



one the tunnel, compare [12|]. This and all other transfer matrices we take as 
a function of the dimensionless variable /3 = Kb, related to the energy and to 
the length of the cell S. We also consider them as a function of the variable 
7 = ub related to the strength of the delta-potentials, and study the family 
of systems with varying strength. Instead of 7 we shall often use the ratio 
5 := % because then the matrix elements of Mi are (linear) polynomials in 
the variable 5. Since we shall generate all other transfer matrices by matrix 
multiplication, a general transfer matrix will be a polynomial in the parameter 
5. This will the basis for a polynomial method for handling the analysis. 
We introduce the following short-hand notation for the elements of a matrix 
Mi 



Mi 



a, 



bi 

d; 



(8) 
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. 1 j 



(<k + di), 
(a, - di). 



The second row of eq. || determines the half-trace %i = \tr{M). For 
matrix M\ we find from eq. [7] 

Vi = -^(Ai-A^ + ^ + Ar 1 ), 
d 1 = \i(l-±S). 
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Figs 2,3: The strings of length Kb — 1, r respectively produces two different 
band germs which encapsulate the same bound state. Here and in all the 
corresponding Figures, the values of 7 and j3 determine the potential strength 
and the energy respectively and are plotted as the horizontal and vertical 
coordinates respectively. 

We shall assume that we are in the band germ and, moreover, that the upper 
edge of the band occurs at negative energy. In sections 4, 5 we shall deal 
with more general scenarios. For < 1, the eigenvalues of Mi are complex 
conjugate of modulus 1 given by 

0\ = x\ + i\Jl — xf, 8 2 = x\ — i^Jl — x\. (10) 

So Xi is the half-trace, and the local Bloch label K for fixed energy is obtained 
from 

cos(Kb) = 7), < K < (11) 

The values (3 : Xi(/3) = —1, 1 determine the two band edges respectively. 
The lower diagonal element by d\{0) = determines, compare |12| and 
what was explained above, the single bound state of the delta-potential at 
5 = ^ Q = 2. The variable Ai, related to the tunnel length, drops out of this 
equation. 

As M is real, the eigenvectors can be chosen complex conjugate. We find 

M X V = VA, (12) 
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with 



v 
P 



V = P I 

V V 

A = * ° 

6 2 

9 1 -a l -yx + idl- x\ 



r 1 f(l - xf) - h(l + Xf) + TiX^y/T^ 



(13) 



With the help of the matrix V we pass to the Bloch system of solutions 
defined by M(x)RV since it obeys M{0)RV = RV, M(b)RV = RVA. In 
line with eq. [12] we mark derivatives with respect to x with a prime and put 

and obtain for < x < b from eqs. |12| , |13"1 , |14| the Bloch state 



Mi{x)RV 



(14) 



P 



>2k 



exp(— KX) 



r 1 



(l-A r 2 )--5(l + Ar 2 ) + 2zAr 1 v / r 



— x\ 



exp{Kx) 



To normalize the current density to the value |^ it suffices to choose det(V) 
pv — vp = —i which together with eq. [13| yields 

<JAi 



(15) 



pp 



P 



1 

%— 

2' 



6X 1 



(16) 



where we have chosen a particular phase for p so that p = —p. The two 
Bloch states obey the bilinear relation 

(-t)(¥ l $i-w l ^) = (-iy + %. (it) 

which can be used to define a scalar product for Bloch states: 
4Def: To any pair $i, <3> 2 of complex solutions of the Schrodinger equation 
we can associate an indefinite hermitian scalar product (<3>i, $ 2 ) by the left- 
hand side of eq. |i~7 . 

The relation between the Bloch states and the exponential states Mi(x)R is 
given from eqs. [L2| and [14] by 



Mi(x)£ = {M 1 {x)RV) V~ l 



(18) 



" $1 


$ 2 




v —p 




$ 2 _ 




—v p 
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A particular case arises if we choose the bound state energy E = — |-k 



and construct the Bloch states eq. 14 for k = k , 5 = 2. This bound state 



2m a 



must increase exponentially in the unit cell towards the delta-potential and 
is given from eqs. [15], [16] and |l~8| for < x < b by 



9 A 

" 1 / $ 1 (x) + $ 2 (x)). (19) 



which differs from the expression given in (12| eq.(22) in the phase for the 
Bloch states and due to the new sequence of tunnel and potential. A natural 
real unbound companion of this bound state from eq. [18] is 

ipi(x) = i($i(x)v - $2(x)v), (20) 

because it can be easily be shown that the scalar product eq. |D] for these 
two real states yields the Wronskian 

= (21) 

Here we wish to comment on the method of orthogonalized plane waves 
OPW in which Bloch states are required to be orthogonal to certain bound 
states [§] p. 206-208. We doubt that this requirement can be handled with the 
standard integral scalar product in Hilbert space, or with a scalar product 
that involves only integration over the unit cell: The first integration does not 
apply without qualification to unbound states, and the second one introduces 
terms at the boundary of the unit cell. 

2.3 The string S n . 

Now we pass from the string S by n-fold repetition to the string S n with the 
transfer matrix Mf . This power of Mi can be written as 

M? = VA n V-\ (22) 

Explicitly the four matrix elements from eqs. 12], |l^, are 



a n = a(M?) = cos(nKb)+ yi Sm ^^ , (23) 

d n = d(Mf) = cos(nKb) - yi Sm( , n 5? , 

sm(Ao) 

,,-„, 1. r sin( nKb) 
2 sm(Kb) 

r (M n )- 1 x -i: sin ( nKb ) 
° n ~ c ( M i)-"2 A i 6 sin{ Kb) ■ 



with yi taken from Mi in eq. [7|. 
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2.4 Rational Bloch labels. 



Consider the transfer matrix eqs. |23| for the unit cell. For fixed integer 
n > at the points 

nKb = firr, fi = 0,±l,...,±(n-l),±n (24) 

it becomes = (— l) n e. For any such Bloch label it follows that the Bloch 
states have the period 2nb with symmetry group Cm- We can relate these 
states to a chain of n atoms by imposing as boundary conditions the factor 
(— 1)™ = ±1 for even and odd n respectively on the interval nb. Moreover 
the Bloch state for fixed /i on the unit cell transforms according to the rep- 
resentation 

£>" = exp(i//^) (25) 
for the generator of order 2n of the cyclic group C 2n - 

There is degeneracy as we have a real hamiltonian: For opposite signs of Kb 
and \i we have pairs of Bloch states with equal energy, running to the left and 
right respectively. The set of rational states for fixed n yields a spectrum of 
(in part) conjugate eigenstates. The energy increases with The rational 
states are shown for n = 2, 3 in Figs. 5, 6. 

With respect to C 2n , each eigenstate is characterized by a real or by a pair 
of complex conjugate irreducible representations of C 2n . All these rational 
states occur in the band and separate it into n partial bands characterized 
by 

K>0: ^r<K< (/i + 1)7r , fj, = 0,1, ... ,(n — l),n. (26) 
nb nb 

and similar expressions for K < 0. 

2.5 Bound states and clusters of the string S n . 

The bound states of the finite string S n are characterized from eq. |3| by 

d(M?) = 0, (27) 
sin(Kb) 



tan(nKb) 



yi 



The free string forms a cluster with n atoms and hence at most n bound 
states. We now compare these bound states with the states of the (up to a 
sign) periodic string of n atoms: In the interval between two successive zeros 
of tan(ni^S) with nKb = fin, (/i + l)7r, which correspond to the rational K 
labels associated with states of the periodic string, the value of tan(ni^6) 
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goes once to infinity. The right-hand side of eq. |27| is a smooth function of 
P except at the point yi{(3) = 0. The left and right-hand side of eq. ^7| are 
displayed for the case n = 10 as functions of (3 in Fig. 4. 




Fig. 4: Left- and right-hand side of the binding eq. [27] on a string of length 
10 for 7 = 4 function of (3. 

On any interval between consecutive zeros of tan(nKb) there is precisely one 
intersection point of the left and the right-hand expression of eq. [27|. 
In terms of boundary conditions this implies that, in between the energies 
of two consecutive rational values of Kb with periodic boundary conditions, 
there is precisely one bound state of the n-atom cluster. This state, termed a 
decaying state in [l2fl, is a single linear combination of two degenerate Bloch 
states on the interval of length nb and, when the n-atom string is cut out, 
admits exponential decay to the right and to the left. 

5 Prop: Any set of rational Bloch labels eq. ^4] separates the band into n 
partial bands. Each partial band contains a single bound state of the free 
n-atom string. 

The relation between rational Bloch labels and bound states is shown in the 
following Figs. 5,6. 
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Figs. 5,6: The strings S 2 and S 3 . In the supercell analysis, two or three band 
germs are glued together without gaps and overlaps. Each band includes a 
single bound or decaying state. The horizontal coordinate is the potential 
strength 7, the vertical coordinate the energy (5. Notice the relation between 
the onset of additional bound states and the subband edges. 

Since the rational labels K form a dense set of energies, we must also have 
a dense set of bound or decaying state energies of free n-atom strings. Each 
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bound state coincides on a cell of length nb with a linear superposition of 



two Bloch states at the same fixed negative energy similar to ip 2 in eq. [LS. 



Each one has a companion state of a form similar to ipi in e Q- POL and these 



two states have a Wronskian as in eq. |2T 



2.6 Participation number. 

From the relation between n-atom clusters and periodic states we can assign 
within a band a participation number n to the rational Bloch vectors of the 



form K — eq. ^J: This number determines the number n of atoms which 



participate in the Bloch states whose period up to a sign is nb, and at the 
same time assigns n bound or decaying states of the cluster. Small values of 
n indicate states which involve a few atoms and may be called localized with 
respect to the clusters. 

2.7 Supercell interpretation. 

An alternative interpretation of the results given above can be given in terms 
of a supercell scheme: We consider the supercell of length nb formed from 
n centers. The corresponding Brillouin zones have the size 27r/(2bn). The 



positive Bloch labels for the n bands in this scheme are obtained from eq. |26 

as 

0<K(n^) = K-^-<— (28) 

nb nb 

\i = 0,l,...,(n-l),n 

and similar expressions for K(n,/i) < 0. With the new Bloch label K{n,u), 
these superbands are identical with the partial bands for the rational values 
of K. 

6 Prop: The superbands for the supercell of length nb are glued together 
without gaps and overlaps and fill up the single band. Their Bloch states are 
identical to those of the partial bands of the rational scheme. Within each 
superband there is included a single bound state. This bound state coincides 
on the supercell with a bound state of the free n-atom cluster which forms 
the motive of the supercell. The alternating pattern of bound and periodic 
n-atom states illuminates the structure of the band. 

2.8 Large n limit. 

Now we consider the limit n — > oo. The rational If -labels form a dense 
set on the Brillouin zone. The length of the periodic chains and of the 



16 



clusters increases with the participation number. As the rational numbers 
form a dense set on the Brillouin zone, we have pairs of Bloch states and 
bound states near any value of K in the band. By picking rational values 
of K characterized by n, we select a set of periodic states which pairwise 
encapsulate bound states of the finite string S n . 

From the equal spacing of the zeros of tan(nKb) with respect to the band 
label K, and from the inclusion property for bound states there follows 
7 Prop: The density of bound and of Bloch states in the limit n — * oo is a 
constant function of the Bloch label K . In this limit, the density of (bound) 
states per unit energy interval fulfills the well-known relation 

dN dK 

Although this is a smooth function in the limit, its rational approximants 
with large but finite n may look quite irregular as functions of (3. 



3 Finite Quasiperiodic strings at negative en- 
ergy. 

As an example of a quasiperiodic system we shall take the well-known Fi- 
bonacci system: We maintain delta-potentials of equal strength for the atoms 
but admit two different intervals S and L of length b and qb respectively. For 
the proper Fibonacci case we choose q — r — |(1 + v^5), but for comparison 
with the periodic case we also consider the value q — 1. 
The Fibonacci system we define algebraically by the recursive words W m and 
initial data in the alphabet (S,L), 

W m+1 := W m ^W m , W l = S,W 2 = L. (30) 

The word W m contains f m letters, f m ~i letters S, and f m letters L. Here f m 
are the integer Fibonacci numbers defined by f m +i = fm-i + fm, fi = fi = 1- 



3.1 Preview: Energy gauge in Fibonacci strings. 

As an introduction we take from (12| the example of two strings S, L which 
represent the same attractive delta-potentials combined with tunnels of length 
b and rb respectively and then pass to the string SL. The energy gauges for 
S and L both contain the same bound state. The string SL has two bound 
states but may have, depending on the tunnel length, two or one band germ. 
In the first case there exists a region of negative energy where the single 
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strings admit band germs but the combined string SL does not. This case is 
schematically represented in Fig. 7. 

f(S) S f(L) L f(SL) SL 



Fig. 7: Fibonacci strings: The string S to the left and its energy gauge f(S) 
are as in Fig. 1. The string L in the middle has the same bound state within 
the gauge f(L) of smaller width. The string SL to the right has two attractive 
^-potentials and hence the same bound states as S 2 in Fig. 1. In contrast 
to S 2 , the string SL has two separate band germs. The corresponding two 
parts of the energy gauge f(SL) comprise each one bound state but block 
the energy of the single-atom bound state. 



3.2 Substitional systems and their invariants. 

The Fibonacci system to be considered can be seen as a particular case of a 
substitution system. We shall restrict our attention to substitutional systems 
generated by automorphisms of the free group. For these algebraic concepts 



in general we refer to fll3| and follow in detail ||10|| . Consider the free group 
F 2 with generators (alphabet) a\,a 2 . Its elements are all the words in the 
generators and their inverses. An automorphism of F 2 is defined as a map 

p:a 1 ,a2^fi(ai,a 2 ),f2(ai,a 2 ) (31) 

which is algebraically invertible so that ai,a 2 can be expressed in terms of 

fli f2- 

Example: Consider the substitution ai,a 2 — > fi = a 2 , f 2 = a\a 2 . Its inverse 

is given by a\ = f 2 fi , «2 = fi- This substitution describes in abstract 

terms the Fibonacci systems to be discussed below. 

For any automorphism of F 2 we have the general theorem 

8 Prop (Nielsen 1918): For p an automorphism of F 2 , the commutator 

JC(a\,a 2 ) := a\a 2 a{ a 2 obeys 

p{K) = wIC^w- 1 (32) 
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with w e F 2 . 

Example: For the Fibonacci automorphism one finds 



p(/C) = KT 1 (33) 

From this abstract algebraic set-up we can pass to matrix systems by first 
mapping the generators ai, a 2 to two elements of a matrix group like SL(2, R) 
for the transfer matrices and then performing the automorphism p induced on 
the matrix group. For certain classes of such induced automorphisms, one 
can obtain the (half-) traces of the images directly from the (half-) traces 
of the generators. This question is studied in the theory of trace maps for 



which we refer to Peyriere [^6 ]. The Fibonacci automorphism and related 
systems belong to this class HT0| . The Nielsen theorem in connection with 
the induced automorphisms now provides one or more invariants of induced 
automorphisms: Consider the induced commutator /C on the matrix group 
SL(2, R) and take its trace. Under the Fibonacci automorphism the com- 
mutator as a matrix is transformed into its inverse. Since all matrices are 
unimodular, we get on the matrix group 

tr{p{K)) = tr{K- 1 ) = tr()C) (34) 

9 Prop: The trace of the commutator is an invariant under the Fibonacci 
substitution. 

The commutator K(gi, g 2 ) in a matrix group is a measure of the non-commutativity. 



In particular if the two matrices commute, the half-trace eq. |34] must be 1. 
Conversely if K, equals the unit matrix, the elements gi,g 2 commute. We 
shall come back to this property in later sections. 



3.3 Recursive calculation of the transfer matrix. 

To represent the Fibonacci system by delta-potentials and the states by trans- 
fer matrices we define 

M(S) := Mi, M(L) := M 2 . (35) 

We represent the Fibonacci string W m+ i of eq. |30| by the transfer matrix 

M m+1 = M m _ x M m . (36) 

The transfer matrix M 2 has the same analytic form as M\ given in eq. ^ 
with the same variable S and the same energy E but with the replacement 

A x = exp(/3) -> A 2 := exp(q(3). (37) 
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Similar systems have been studied extensively in the literature, compare 
references given by Kohmoto || and in fl3|. Many of these studies used 
a discrete version or examined the recursion numerically. In what follows 
we shall use the recursive method to explore the analytic structure of the 
Fibonacci strings as functions of the variables (3,5. We recall that Mi,M 2 
eqs. |7|, |35] are linear polynomials with respect to 5. From this property and 
from eq. [36| we get immediately : 

10 Prop: The matrix elements of M m are polynomials of degree f m with 
respect to the parameter 5, with coefficients which are functions of /3 via the 
expressions Ai, A 2 . 

For the finite Fibonacci string we wish to study the matrix elements of the 
transfer matrix Mj and in particular their combinations Xi,yi,di = (xi — 
Hi), since they yield information on the eigenvalues and bound states. Two 
methods are available for the computation: 



(i) We can use the full matrix recursion eq. [3(| with the starting matrices 
Mi,M 2 . 

(ii) We can use recursion techniques for the half-traces and the other matrix 
elements, as was proposed in [fl(|. From this reference we deduce the following 
recursion equations for the combinations of matrix elements: 



a m+l b m+ i 



Vm+1 
dm+l 



Q>m d m 

a m -\- d m 

—C m -2 d m -2 
2,X m X m —i %m— 2i 
^•EmDm—Y ~\~ Vm—2i 
•Km+1 Vm+1- 



Q"m—1 bm—l 
Cm— I d m —\ 



(38) 



To start these recursive relations we need the matrices Mi,M 2 given from 
eqs. 0, [35| and compute M 3 from 



M 3 = M X M 2 : 

a 3 = Ar 1 A 2 1 (l + i,5) 2 -iA 1 A 2 1 ( 5 2 , 

h = ^A^l + i*) + ^A^l - \5) 

c 3 = -Ia^a^i + ^-^a^i-^) 

d 3 = XiX2(l-\5) 2 -^ l \ 2 5 2 . 



(39) 
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These matrix elements are polynomials of degree 2 in the variable 5 as ex- 
pected. The half-traces x m form a recursive system by themselves. This is 
the well-known Fibonacci trace map. 



11 Prop: By use of the recursive relations eq. [3^ we can construct algebraic 
polynomial expressions of degree f m for the half-trace x m (/3,8) and for the 
matrix element d m (/3,8), which we call the band polynomial and the bound 
polynomial respectively . The edges of the band germs for the Fibonacci string 
W m are given by the f m roots of the band polynomial equations 



x m (p,5) = ±1, 



(40) 



while the bound states are given as the f m roots of the bound polynomial 
equations 

d m (P,5) = 0. (41) 

The f m roots of the polynomials eq. |I] are real since they represent special 
half-traces of real matrices. The f m roots of the polynomial eq. |T] are real 
for sufficient strength of the delta-potential since then we must get f m bound 
states. 

For low values of m the roots of the polynomials can be obtained in closed 
algebraic form. 
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Figs 8-11: Band germs and bound or decaying states as functions of 7 hor- 
izontal and (3 vertical for Fibonacci strings of 3-6 atoms. The bound states 
are not always inside the band germs. The band germs which were glued in 
the periodic string are now all separated. The number of band germs increases 
proportional to the number of atoms. 

In Figs. 8-11 we use the roots of eqs. HO, PHI to plot for the first 6 Fibonacci 



strings the edges of band germs and the bound states as functions of (3, 
negative vertical coordinate, and 7, horizontal coordinate . For each value of 
m, we compare the periodic case q — 1 with the quasiperiodic case q = r. 
For the value q = 1 we are back at periodic cases and would obtain glued 
systems of 2, 3, 5, 8 superband germs respectively in a way similar to Figs. 
5,6. These superband germs encapsulate the bound states of the free clusters. 
For the Fibonacci case q = r there appears a system of f m band germs, for 
any value of 7 separated by gaps. The bound states for higher values of 7 are 
encapsulated within these band germs. For large fixed values of 7 the band 
germs become narrow and form a triple of subsets. The latter feature may be 
related to the occurrence of isolated 2-atom clusters at the short distance b in 
the Fibonacci strings. These may be considered as examples of motive clus- 
ters discussed in section 1. The locally symmetric or antisymmetric hybrid 
states of these clusters can be related to the lowest and highest subset. 
Figs. 8-11 also demonstrate how the quasiperiodic system develops under 
the Fibonacci recursion or inflation. If one looks for fixed strength 7 at 
the interchange of bands and gaps, one observes a permanent change of the 
density of states (DOS). In an approximant computation of this quasicrystal 
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one would stop at a certain length and obtain the DOS from the band germ 
analysis of the corresponding string. In view of the present examples it is 
hard to believe that this method describes the real DOS of the quasicrystal. 

In a periodic potential we get a finite number of bands at least at negative 
energy which can encompass a number of electron states proportional to the 
number of atoms. In the Fibonacci strings notice that, in marked contrast to 
the periodic case, when extending the Fibonacci string the number of band 
germs which encapsulate bound states increases proportional to the number 
of delta-potentials, i.e. proportional to the number n of atoms, each with one 
electron, in the model. If this is so, then the DOS cannot be computed from 
the gap and band germ structure of the string. As we have only n electrons, 
it makes no sense to distribute them into a system of n bands. We conclude 
that in approximant computations of the DOS one should keep track of the 
number of bands in relation to the system size. Otherwise the conclusions 
on the DOS in the proper quasiperiodic system are doubtful. 

Now we turn to the significance of the conserved quantity related to the trace 
of the commutator. This trace is related to the invariant I considered in || 
and i by 

-tr{K) = 21 + 1. (42) 



If the commutator becomes the 2x2 unit matrix we get 1 = 0. In |12| we 
computed the commutator for the present Fibonacci system with the result 

1 1 11 

-tr(/C) = l + -(-) 2 (sinh( K (r-l)6) 2 (43) 

2 Ik 

This expression is always larger than 1. It implies that in the negative energy 
tight binding scenario the two transfer matrices never commute. A different 
result will appear for positive energy and scattering. 



4 Periodic strings at positive energy. 

We turn to electron states at positive energy. To this end we use the analytic 
continuation from the variable k to the variable k defined by 

K — ► —ik. (44) 

For the energy we have now E = ^k 2 > 0. For the variables /3, 5 we shall 
use 

p -> -i/3, 5 -> iS. (45) 
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So we maintain the symbol (3 as a real variable in what follows. We shall 
plot the positive energy range E = ^j^P 2 at negative values of the new real 
parameter (3. Inserting the analytic continuation into the transfer matrix M\ 
eq. we find 

x x = cos(p) - ^6wa.(J3), (46) 

yi = isin(/3) +i^5 cos(fi), 
Ai = exp(— i/3) = exp(-ikb). 

Similar expressions apply for M 2 - 



4.1 The <S-matrix. 



The transfer matrix of a finite string can be rationally related to the scat- 
tering matrix S. For this purpose we shall use the transfer matrix in a form 
similar to eq. [5] which corresponds to exponential functions. The analytic 
continuation eqs. |4"4] , |4"5| is made both in M(x) and with the matrix R := -^R 

used instead of R. We redefine for positive energy 



M(x) := R~ 1 M(x)R. 



(47) 



The negative energy tunnels with M given by eq. [| become positive energy 
channels. The new fundamental system of solutions is 




exp(ikx) exp(-ikx) 
ikexp(ikx) —ikexp(-ikx) 



(48) 



and consists of free plane waves running to the left and right respectively. 
The transfer matrix in the channels becomes 



M{x) := R~ l M(x)R 

exp(ikx) 
exp(— ikx) 



(49) 



Consider now a finite string of length h with transfer matrix M and with free 
channels on the left and right respectively. Denote the amplitudes on the left 
and right of the string by l+,l- and r + ,r_ respectively. These amplitudes 
are related by 

», i r ^ i, i r ; 

(50) 



r+ 




a b 






r_ 




c d 
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The elements of the scattering matrix are now determined by the ratio of 
amplitudes under certain boundary conditions and by phase factors which 
account for the length h of the string. For scattering from the left we obtain 



cl + + dl- = 0, 

— - exp(-ikh) 
H- 

d~ l exp(-ikti), 
Z_ 

h 

-d- l c. 



(51) 



For scattering from the right we get 



l + = dr + — 6r_ = 0, 

S = — exp(-ikh) 

r_ 

= — gT 1 exp(— ikh), 
S + - = — exp(— 2ikh) 



d bexp(—2ikh). 



(52) 



Here S ++ and 5*_ + are the amplitudes of forward and of backward scattering 
respectively. These expressions are given and discussed in [[[(]] eqs. (50-54). 
12 Prop: The scattering matrix S has elements which are rational expres- 
sions in the transfer matrix M. The full scattering matrix becomes 



5_ 



Si 

5_ 



(53) 



d 1 exp(-ikh) d 1 bexp(—2ikh) 
—dr x c d^ 1 exp(-ikh) 



The exponential factors arise by requiring that for the free transfer matrix on 
the string of length h we must have S = 1. The expression for the iS-matrix 
in terms of the transfer matrix is non-linear but rational. We must compute 
S from the transfer matrix of the full string, or otherwise evaluate all the 
backscattering effects from the components of the string. The unitarity of the 
iS-matrix and its transformation under time reversal follow from properties 
of the transfer matrix ||10|| . 
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Consider now the analytic continuation eq. 



in the inverse form 



k — ► %K, K > 0. 



(54) 



From inspection of eq. [53| under this analytic continuation we find 

13 Prop: The poles of the scattering matrix eq. [53| on the positive imaginary 

/c-axis are given by cI(k) = 0. These poles determine the bound states of the 



system, fully in line with the analysis given in eq. 41 



4.2 The <S-matrix for the periodic string S n . 

Again we consider a range of positive energy such that \x\\ < 1. This range 
may be called a band germ although it cannot correspond to bound states of 
the system. In this range we can find the eigenvalues and eigenstates of M\ 
by analytic continuation of the expressions eqs. ||, || . For the periodic string 
S n with the transfer matrix M™ we find the same expressions as in eq. |3|, 



but with yi and Ai now taken from M\ in eq. |46|. We shall again introduce 
a band label K by 

cos(Kb) = xi = cos(P) — -<5sin(/3). (55) 

Consider first the backward scattering determined from eq. [53| by c n and d n . 
The matrix element c n oscillates rapidly with large n. We obtain maximum 
backscattering at the discrete values of the -fT-label 



-i-Ar^-l)"^"" 1 ^ 



Kb = //7r,// = 0,±l, (56) 

d n = (-l)^(l_ nyi (_l)M) ) 

exp(-ikh) = exp(-inkb). 

From 

cos(//7r) = (-1) M , (57) 

these maxima correspond to the edges of the positive energy band germs and 
become sharper with increasing n. 



A high-energy limit is obtained for j3 = kb 3> 1. With 5 = 7//3, eq. |55 
enforces 

p = kb = Kb + z/2vr + e. (58) 



Together with eq. [56|, this equation yields discrete values of k on the points 



of the reciprocal lattice. 
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14 Prop: For large n and sufficiently high energy, the backscattering ampli- 
tude takes its non- vanishing values on the points of the reciprocal lattice. 
To lowest order in e we find 



*i = (-1)", (59) 

d n = (-ir(l -i~n6). 
In the high-energy limit we then get from eq. |53] at the band edges 

s ++ = ids? (60) 



S. 



'2 

i\n8 
1 — ihn5 



The limits n — ► oo and of large energy should be distinguished from one 
another. In the limit n — > oo the absolute value of the backward scattering 
amplitude approaches its maximum \S-+\ = 1 allowed by unitarity, while the 
forward scattering amplitude goes to zero. 



5 Quasiper iodic strings at positive energy. 

The analytic continuation of the transfer matrix works for the Fibonacci 
strings. We simply must start the string with the analytic continuation of 
the starting matrices M 1? M 2 as given for M\ in eq.09. Then we can apply the 



recursion technique of eq. [38] to find an algebraic expression for the transfer 



matrix of the Fibonacci string W„ 



5.1 The Fibonacci-atlas. 



In principle we can construct the elements of the transfer matrix as poly- 
nomials of degree f m with respect to 6 for any order m. We obtain the 



corresponding iS-matrix from the relations eq. [53|. In general it will be hard 
to obtain in this way closed limiting expressions as was possible in the peri- 
odic case. Nevertheless there are particular regions where closed expressions 
can be derived. We illustrate these regions in Fig. 12 by plotting the band 
germs of both transfer matrices M 1; M 2 for q = r as functions of the variables 
/3,7 in the range (3 > 0, 7 > 0, 7 < 0. This plot we call the Fibonacci- atlas. 
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Fig. 12: Fibonacci atlas: Band germs of the transfer matrices Mi,M 2 for 
positive and negative energy as functions of the energy variable (3 and the 
strength variable 7. Upper band edges drawn as full lines, lower band edges 
as dashed lines, heavy lines for the matrix M 2 . The two weak horizontal 
dotted lines mark values (3(p) where the two transfer matrices commute and 
give rise to extended states. 

Consider again the commutator K of the two transfer matrices. We must 



insert the analytic continuation to positive energy and obtain from eq. 43 



1 1 ?/ 

-tr(K) = l + -(^n S m(k(r~l)bf 



(61) 



This expressions shows that there are periodic points where the half-trace 
becomes 1 and the invariant becomes 1 = 0. These points are given by 



(3{p) = rpn, p = 0, 1, 



(62) 



Moreover it was shown in § that the deviations I — 1 at these points starts 
quadratically with (/3 — (3{p)). The values (3{p) are marked by dotted horizon- 
tal lines in Fig. 12. If a line of this type intersects, depending on 7, with over- 
lapping band germs for both transfer matrices, we obtain from A 2 = (— l) p Ai 
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their proportionality, 

M 2 = (-1) P M X . (63) 



15 Prop: At the values eq. |6T], the two transfer matrices M 1; M 2 commute 
and are even proportional to one another. 

More precisely we require three conditions at these special points: The half- 
traces of both transfer matrices should be smaller than 1, and the commuta- 
tor should become the unit matrix. These conditions are controlled by the 
Fibonacci atlas Fig. 12: Consider a dotted horizontal line corresponding to 
commuting transfer matrices. For not too large values of 7, it runs inside 
overlaps of the two band germs which imply that both \x\\ < 1, |x 2 | < 1. 
Outside this region, that is outside the points marked by circles in Fig. 12, 
we still have commuting transfer matrices, but their individual half-traces 
are larger than 1. It follows from algebraic properties of the group SU(1, 1) 
given in [|l(J] that a vanishing commutator implies group elements of the same 
class type. 

We can now see the advantage of the present algebraic and polynomial ap- 
proach compared to numerical studies of the trace systems: For example in 
Q the points eq. ^2] were explored numerically and for fixed values of a 
strength 7 < 0, corresponding to repulsive delta-potentials. The Fibonacci- 
atlas of Fig. 12 now yields from closed algebraic expressions a full view on 
the regions of commutativity for all positive and negative values of 7 and 
positive energy. 

In a neighbourhood of such a commutative point we get, using the commu- 
tativity, for the full Fibonacci string W m the transfer matrix 

M m = (-iy^M[ m . (64) 

This transfer matrix agrees up to a factor with the transfer matrix of a 
periodic string of f m delta-potentials with the spacing b. It follows that the 
iS-matrix of the Fibonacci strings near these points has the form and limiting 
values discussed in section 4.2. 

16 Prop: In the neighbourhood of the discrete positive energies correspond- 
ing to eq. |62|, the transfer and 5-matrix of the Fibonacci system are equiva- 
lent to those of a periodic string of length f m b. 

Some care is required because /3(p) in eq. ^ is not rational whereas for the 
periodic case we used the discrete labels Kb and j3 = kb. Note however that 
the value /3{p) for p = fi and not too small values of / is well approximated 
by the integral multiple (3 = //+i7r of tt. 

Finally in Fig. 13 we give the real wave function of an electron travelling at 
positive energy through a Fibonacci string of attractive delta potentials. The 
energy is tuned to a commutative value. At each passage through a delta 
potential, the derivative of the wave function jumps by a finite value. 
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Fig. 13: Electron state propagating through a Fibonacci string of attractive 
delta-potentials at a positive energy where the transfer matrices commute. 

6 Conclusion. 

Our study of a very simple ID electron system has demonstrated most of 
the points stated in general terms in the introduction, section 1. The local 
boundary conditions were implemented and related. The negative and pos- 
itive energy scenarios showed different behaviour. At negative energy, the 
bound states have a clear relation to the Bloch states. In going from peri- 
odic to quasiperiodic strings, the glued and gap-less systems of superbands 
open gaps and split into subbands which encapsulate bound states. Motive 
clusters appear in the bound state energy. With increasing length of the Fi- 
bonacci string, the subband structure changes in a non-trivial fashion, with 
a number of subbands proportional to the number of atoms, and so puts 
question marks to approximant calculations. At positive energy, the band 
edges in periodic strings are related to maxima of scattering amplitudes. For 
quasiperiodic Fibonacci strings we gave closed algebraic expressions. At pos- 
itive energies related by a period wrt. k, the transfer matrices commute and 
give rise to maxima in the scattering amplitudes. 
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